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Abstract 

It is known that an accelerating charge radiates according to Larmor formula. On the other 
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qj^ hand, any DC current following a curvilinear path, e.g. a circular loop, consists of accelerating 
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charges, but in such case the radiated power is 0. The scope of this paper is to analyze and quantify 
how the radiation vanishes when one goes to the continuum DC limit. 
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I. INTRODUCTION 



There are many physical configurations in which discreet charges are in acceleration, 
but in spite of that they behave almost like steady state DC. This situation is encountered 
in any DC or low frequency electrical circuit, because in spite of the fact that the charge 
distribution is always discreet and each charge accelerates in the influence of the electric 
field, one may consider the charge distribution as almost continuous. For example the 
model used for conducting materials is of an average drift velocity for positive and negative 
charges which is proportional to the electric field v 1 * 1 = ±/i ± E, /x ± being the mobility of the 
positive and negative charges, respectively. Hence the current density is J = p + v + + p~v~ = 
(p + p + — p _ /i _ )E, where p is the charge density of the positive and negative charge carriers, 
p~ being negative, so that the total charge density is p = p + + p~ = e V ■ E. In case the 
charges are protons and electrons, in a conductor, the protons mobility is p + = and 
p + = —p~, so that usually V ■ E = like in free space. In addition, the model considers p _ 
to be almost uniform, so that one defines the conductivity a = —p~p~, resulting in Ohm's 
law J = o~E. 

Another situation, for which charges are somehow "more discreet", are DC or low fre- 
quency ion drift devices, and for those we usually have one type of charge carriers, say 
positive ions. Here we use J = pv = ppE, and one may not assume the charge density 
is uniform, but rather has to use Gauss's law eoV • E = p, resulting in a set of nonlinear 
equations. In principle, such problem is time dependent, because discreet ions move in the 
space, but it appears that the DC approximation V • J = works very well for such cases 
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In either of the above situations, currents may follow curvilinear paths, in which case, the 
charges clearly accelerate, also if the magnitude of their velocity is constant, and still the 
DC model works well for most practical cases for which the density of the charge carriers is 
big. 

The purpose of this work is to understand how the radiation vanishes when the density 
of the charge carriers approaches the continuum steady state. Certainly, the radiation has 
io disappear gradually, so that this vanishing may be quantified. Some preliminary work 



10( has been done in this direction, but because this work has been done a priori in a non 



relativistic approach, its results are inaccurate and incomplete. 
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The calculation is done in a canonical configuration of charges in circular motion at con- 
stant speed. The configuration and the formulation are explained in section 2. In section 3 
we calculate the fields, explain their behavior and derive an exact expression for the radi- 
ated power. For completeness, we also calculate the radiation reaction and show that the 
power needed to support the radiation equals the radiated power. In section 4 we derive an 
asymptotic result for the case the number of charges is big (continuum limit). As a special 
case, we also derive the limit for slow charges, and this case is of importance because it 
represents the typical situation of DC currents in devices as discussed before. The paper is 
ended with some concluding remarks. 

The work is written in SI units, and we shall use the known constants: vacuum permittiv- 
ity e = 8.85 x 10~ 12 F/m, vacuum permeability fiQ = An x 10~ 7 H/m, free space impedance 

= a/^oAo = 376.73 f2 and speed of light in vacuum c = 1 / ^Uo e o ~ 3 x 10 8 m/sec. 

II. CONFIGURATION AND FORMULATION 

A total amount of charge q is rotating in a circle of radius d at constant speed v so that the 
angular velocity is u — v/d. The charge q is "split" into N charges of value q/N, uniformly 
distributed around the circle, so that the charge number k is at the angle ut + 27rk/N, where 
k = 0,1,.. A - 1. 

The configuration is shown in Figured], for N = 3. 

The location of the charge k as function of time is given by: 



The fields propagate with the speed of light c. Hence the fields at the observer location 
r at time t are influenced by the motion of each charge, at an earlier (retarded) time. 
Specifically, the fields are influenced by the motion of the charge k at time t' k so that 



r' k (t) = d[xcos(tut + 2nk/N) + y sin(wt + 2nk/N)\ 



(1) 




(2) 



At large distance from the charges, one may approximate: 



R k pa r — d sin 6 cos 



(3) 
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(r,8,cp,t) ^ 




FIG. 1: (color online) iV charges (here N = 3) of magnitude q/N each, in circular motion on the 
xy plane at radius d around the z axis. The far observer's spherical coordinates are also shown in 
the figure. 



where 



<f) k = ut r k + 2-nk/N - if, 
hence the retarded time t' k may be calculated from the following implicit equation: 



(4) 



t' k = t — r/c+ (d/c) sm9 cos4>k, 



(5) 



which may be solved numerically by setting a "1st guess" t' k = t — r/c in the right side 
of eq. (JSJ) and recalculate t' k until convergence is obtained. 

Figure |2] emphasizes the meaning of retarded positions of the charges. 

Our purpose is to calculate the power radiated by those rotating charges. Of course, for 



N — 1 the result is known by Larmor formu 
The general form of the Larmor formula 



a (and will be confirmed later on). 
J-lal is: 



Larmor 



g 2 7 6 
67re c 



2 - 09 x /3) ; 



(6) 
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Far observer at 



Charge 1 at time t 
Charge 2 at time t 




(r,8,cp,t 



Charge 1 at time t\ 



■■'&: 



^ 

Charge 2 at time t 



Charge at time t 
^ Charge at time t' 



FIG. 2: (color online) The dark colored spheres represent the charges at the current position at 
time t and the light colored spheres represent the charges at the retarded positions at times t' k . The 
retarded positions are connected with dashed lines to the observer location, and those distances 
are called Rq, R\ and i?2- This emphasizes that the field at observer is determined by the velocities 
and accelerations of the charges at the retarded times. 



where 



(3 = v/c. 



(7) 



and 



7=l/Vl-/3 



(8) 



If one defines the angle between the velocity and the acceleration as a, one may express 
((3 x /3) 2 = (/3/Ssina) 2 an rewrite 



PLa rm or = f^/3 2 (l-/3 2 sm 2 a). (9) 

In our case of circular motion, the velocity is perpendicular to the acceleration, so that 
sin 2 a = 1. Therefore by replacing 1 — (3 2 = I/7 2 , the Larmor formula simplifies for our case 
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to 

PI nm 

ljv=1 " 6ne c ~ 67re oC 3' [W) 

where a is the acceleration. We will calculate how this power decreases when the number 
of charges N increases. 



III. FIELDS AND POWER CALCULATION 

To calculate the power radiated from the collection of charges in Figure HJ one needs only 
the far fields, i.e. those who behave like 1/R. The far electric and magnetic fields due to 
the moving charge k are given by 

4tt(1 - (3 k ■ R k yR k 

and 

H fc = K k x E k / Vo (12) 

where R& is the unit vector pointing from the position of the charge to observer and R k 
is the distance between the charge and the observer, as defined in eq. ([2]) - see Figure [21 
/3 k = Vfc/c = r' k /c is the velocity relative to c and a fc = v fc is the acceleration of the charge. 
All the dynamical variables are evaluated at the retarded time (defined in eq. (jSJ)). 

Defining r as the unit vector pointing from the coordinates origin to the observer, we 
may calculate in the far field the difference: 

R*-?= 1 ^4t--«-^ (13) 
| r — r k | r r 

Hence one may use r instead of in eqs. ( ITT1) and ( Tl2l) with an error of order 1/R 2 , 
which does not affect the calculations of the radiated power. Also, in the denominator of 
eq. (ITT]) we may set Rk = r, as always done for far field. So we express the electric and 
magnetic fields as the sum of the contribution from all the charges: 
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AT— 1 



E _ m r x [(r - f3 k ) x a k ] 

AnrN^ (1-/V?) 3 1 ) 



and 



H = r x E/r^o (15) 
Now using eq. (Op), we evaluate r x [(r — (3 k ) x a^] and express it in spherical coordinates: 

r x [(r — (3 k ) x a k = a[6 cos 6 cos (f>k + (p(/3sm6 + sin0fc)] (16) 
and f3 k ■ r evaluates to 

(3 k ■ r = — (3 sm.6 sm.(f) k . (17) 
Setting those results into eq. f[T4l) . we obtain 



E= i^-[£cos0Fc + £Fs] (18) 



and 



H = qa [-0Fs + £cosflFcl (19) 



where the functions Fc and Fs are defined as: 



iV-l 



Fc(t, V ,6,P,N) = J2 ic (<Pk) (20) 

fc=0 

and 

iV-1 

Fs(t,ifi,e^,N) = J2H<l>k) (21) 

fc=0 

and the functions fc and fs are defined as 



fc (0fc) = 77"7 ■ 7 x V ( 22 ) 



(1 + psin 

and 



fs(0fc) = 



p + sin <j) k 
(1 +psin fc ) 3 ' 



(23) 



and p is defined as 



p = (3sm9 



(24) 



and is a parameter which controls the behavior of <fik, as will be soon shown. 
The power per unit of normal area (or Poynting vector) is given by, S = E x H = rE 2 /rjo 
which results in 



S = r 



qa 



AircNr 

The total power is calculated via 



Vo 



(p Fs + 6 cos 9 Fc 



(25) 



2tt pit 

P = r z I dip d9sin9S-r 
'o Jo 



which results in 



(26) 



P=P\ n=i Git, (3, N) = p-^G(t, 0, N). (27) 

Here we factored out the Larmor formula for the radiation of a single charge - see eq. ffTOl) . 
so that G(t, (3, N) is dimensionless and represents the decay of the power. The function 
G(t,(3,N) is given by 



G{t,f3,N) 



277 



f'TT 

dcp / d9 sin9 F(t,<p,6,p,N) 
Jo 



87r 7 4jV 2 J 

hence G = 1 for N = 1, for any t or (3 and the function F is 



(28) 



F(t,ip,9,(3,N) 



ipFs + 6cos9 Fc 



Fs 7 + cos 2 9 Fc 2 



(29) 



Now we eliminate t' k from eq. ([51) and rewrite the implicit eqs. (JSJ) and (J4j) in terms of <pk 



(pk = u(t — r/c) — ip + 2irk/N + p cos <pk 



(30) 



This allows us to change variable ip' = tp — oj(t — r / c) in ( 1281) obtaining: 



8 



G(t,/3,N) = ^ / ( /} dip' [ ^sinfl(Fs 2 + cos 2 flFc 2 ), (31) 
so that eq. ( 1301) is rewritten as 

fc = —if' + 2nk/N + p cos0 fe , (32) 

We see that if and ip' satisfy eq. (1321) . also <pk ~ 27r and ip' + 2tt satisfy it, hence COS0& 
and sin <pk are periodic functions of ip', with a periodicity of 2ir. Therefore, the dip' integral 
in eq. fl3T|) may be evaluated over any period of 2rr, showing that G (and therefore also the 
radiated power P) does not depend on time, so that we may simplify eq. (l3Tj) to 

G(0,N) = - — -n^r / dip' d6sm6 (Fs 2 + cos 2 #Fc 2 ). (33) 

We redefined G to be time independent, so that fc, fs, Fc and Fs in eqs. (I22|) . (123~1) . (!20|) 
and (l2Tj) become functions of instead of i and <£>. 

For understanding the behavior of Fs and Fc, we plot fc, fs in eqs. (122"]) and ( 123|) . and 
theirs sums (eqs. fTSU]) and f|2T]) ) for different parameters. 

Clearly for very small p in eq. (I3"2"j) . 0^, w — <// + 2nk/N, hence the cosine or sine of 4>k 
equal approximately to the cosine or sine of — ip' + 27rk/N, so that both have harmonic 
shapes as function of ip'. In such case fc ~ cos(— ip' + 2irk/N) and fs ~ sin(— ip' + 2-irk/N), 
hence they sum to a small value as observed in Figures [3] and HI For N = 3 (Figure [3]) 
the amplitudes of Fc and Fs are around 0.1, and this decreases with N, as may be seen in 
Figure IU for iV = 10. 

As p increases, fc and fs in eqs. ( 12"2"|) and fl2"3"|) get more distorted, hence they sum to 
bigger amplitudes, as observed in Figures El [61 El and EJ 

In Figure [51 showing the behavior for p = 0.3 and iV = 3, Fc and Fs have amplitudes of 
0.87 and 0.76, respectively, and those decrease for iV = 10 (Figure [6]) to 0.0086 and 0.0082. 
Also we see that for p = 0.3, although fc and fs are distorted, the sums Fc and Fs are almost 
undistorted, unlike for the p = 0.5 and N = 3 case in Figure [71 

The case of p = 0.5 is shown in Figures [7] and [U In this case, not only fc and fs are 
distorted, but also the sums Fc and Fs, however the distortion of the sum decreases when 
the number of charges iV increases, as may be seen for the case N = 10 in Figure El In the 
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FIG. 6: (color online) The fc in eq. (|22p and their sum is shown in panel (a) and the fs in eq. (|23l) 
and their sum is shown in panel (b) for N = 10 and p = 0.3. 

last case the amplitudes of Fc and Fs are 0.59 and 0.5, much bigger than in the parallel case 
with p = 0.3 . 
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FIG. 7: (color online) The fc in eq. (|22p and their sum is shown in panel (a) and the fs in eq. ([23 
and their sum is shown in panel (b) for N = 3 and p = 0.5. 



To summarize, the functions Fc and Fs increase with p and decrease with N, tending to 
undistorted harmonic functions, for big values of N. 

It is interesting to remark that the average of the functions fc and fs is always 0, although 
this is not always visible for the fs functions. This may be proved by calculating 



(fc,s) 



1 

2^ 



2tt 



d<//fc,s( 



1 

2^ 



2- 



rf /_ hc,s(0 fc ) 



(34) 



'o Jo l 1 +psm0 fc ) 3 

where for brevity we called fc,s the functions fc or fs, and we called their average (fc,s) 
We also use the abbreviation hc,s for the functions he and hs denned as 
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FIG. 8: (color online) The fc in eq. (|22p and their sum is shown in panel (a) and the fs in eq. (|23l) 
and their sum is shown in panel (b) for N = 10 and p = 0.5. 



hcOfc) = cos0 fc 



(35) 



hs(0fc) = p + sin. 



(36) 



for the fc and fs average calculation, respectively. We change variable from ip' to <pk, and 
we find from eq. (132]) that 



dip'/d<f) k = -1-p sin^jfc, 



(37) 



getting: 



(fc,s) = - 



2^ 



0fc(O)-27r 



^fc(O) 



hc,s( 



'1 +psin. 



(38) 



where 0^(0) is the value of 0^ at ip' = 0. Because the integrand has a periodicity of 2n, 
one may integrate over any period of 2tt, obtaining: 



1 



(fc,s) = - 

Z7T 



2- 



hc,s(0 fc ) 



'1 +psin0 fc ) 2 



(39) 



This integral may be solved by the residue method on the complex plane. After changing 
variable z = exp(i0fc), one obtains 
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(fc,s) 



1 2i 
2tt p 2 



dz- 



lc,s(zl 



1 2i 



c \z - z x ) 2 {z - z 2 ) 2 2np 2 Jc 



dz Lc,s(z) 



where C is the counterclockwise unit circle integration contour shown in Figure 



(40) 
and 



lc,s are abbreviations for 



lc(z) = z 2 + 1 



and 



(41) 



ls(z) = -i{z 2 - Abz - I) 
where b is a pure imaginary number defined by 



(42) 



b = —0.5ip, 



and the 2nd order poles are (expressed in terms of p or b) 



Zl,2 



-i (l/p =F ^(1/p) 2 - l) = i (-6- 1 ± v/F^+4 



(43) 



(44) 



where indices 1 and 2 refer to upper and lower signs respectively - see Figure fl9]). 
In the last expression of the poles in terms of b, the magnitude under the square root is 
real and negative, so the square root is understood to be positive pure imaginary. 
Also we named the integrand in eq. (j4"0|) Lc,s, defined as 



Lc,s 



\c,s(z) 



[z - z{) 2 (z - z 2 ) 2 



(45) 



to ease manipulations. We remark that only z\ is inside the integration contour (see 
Figure (J9j)), and its residue is 



Res(L, zl) = lim — [Lc,s(^)(2 — Zi) 2 ] = — 
z^z! dz dz 



\c,s(z) 

(z - z 2 y 



Z = Z\ 



h,s (z 1 )(z 1 - z 2 ) - 21c,s'(zi) 
(zx - z 2 f 



(46) 

Using the relations z\z 2 = —1 and z\ + z 2 = —2i/p, one finds that for both fc and fs cases 
the nominator of eq. (|46p is 0, and hence 
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Almag{z} 




Real{z} 



FIG. 9: (color online) The complex z plane on which we show: the poles Z\ 2 in eq (|44p . the 
integration contour C used in eqs. (|3D|) and (f96|) . and the integration contour Ci used in eq. f 1 1 1 [) . 
The poles Z\2 are negative pure imaginary, and \z\\ < 1 and \z%\ > 1, the integration contour C is 
on the unit circle and the integration contour C\ is on a circle of radius smaller than \z\ \ . 



<fc> = (fs) = 0. 



(47) 



Now we continue with the calculation of G in eq. Rewriting eq. for the charge 

m instead of k results in 



—if' + 2-Km/N + p cos ( 



(48) 



which may be rewritten as 



(f) m = -(<£>' - 2n(m - k)/N) + 2nk/N + p cos< 



showing that if we could explicitly express (</?') from eq. 
function of only shifted: 



(49) 



would be the same 
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M<S) = M<p' - M™ - k)/N) (50) 

as also evident from Figures [3] -[HI Therefore both Fs and Fc functions remain unchanged 
for a shift of multiples of 2tt/N, say: 



- / / , 0nr/M \ n a aa cos0 fc (y/ + 27r/AQ COS0 fc+1 (^) 

and the sum being on all charges (and k is modulo AT), we are left with the same result. 
We may therefore rewrite in eq. fl33|) 



/ v = E / ^' ( 52 ) 

JO n=0 J2nn/N 

and after changing variable y?" — tp' — 27m/N, we are left with N identical integrals over 
the period to 2ix/N . For simplicity we rename <p" back to tp' and rewrite the function G 
as: 

o p27r/N PIT 

G(8,N) = j— / dtp' d9sm9 (Fs 2 + cos 2 # Fc 2 ), (53) 

87T7 4 iV y J 

and this will significantly reduce the time of a numerical integration. Now looking at the 
9 dependence of Fs and Fc, we remark from eq. (132]) that <pk depends on p = (3 sin 9, hence 
it is invariant under replacing 9 by 7r — 9, and so is cos 2 9. We may therefore replace the 
integration from to tt by twice the integration from to tt/2, getting 

G(/3,N) = — / dtp' d9sin9 (Fs 2 + cos 2 #Fc 2 ), (54) 

47r7 4 7V J J 

Now we change to the variable p defined in eq. ( )24l) and rewrite eq. ( 1541 obtaining: 



G(A ^ ) = 4dwl d(P 'J dpP ( FS2/ ^ 1 " WW + Fc Vl - (P//?) 2 ) , (55) 
For the case of N = 1 we know G must be 1, but for consistency we shall prove it: 
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G(f3, 1) 



4?T7 4 /3 2 



2tt 



dip' / dpp 
ii Jo 



p + sin O 
(1 + psin (f) 



cos </> 
(1 +psin 0o ) 



(56) 



where Fc 2 and Fs 2 reduced here to a single term. By changing integration order, we may 
perform the dip' integration by the change of variable dip' jd^^ defined in eq. (I3"T|) . obtaining 



G(p,l) 



4tt7 4 /3 2 



(p + sin 0o ) 



^0o(O) 
dpp / G?0o . , 

n ^ (o)-2tt V(l +psin0 o ) 5 v /l - (p//3) 2 



cos 2 ov / l-(p//3) s 



(1 + psin( 



(57) 



All the functions have a 2n periodicity on O , so one may use any limits of interval 2tt for 
0o- By changing variable z = exp(i0 o ) and using the residue method on the complex plane 
we obtain 



2~ 



COS 0o 



ii 4 + p 2 



(l+psin0 o ) 5 4(l-p 2 ) 7 / 2 



(58) 



and 



2rr 



(p + sin0 o ) 2 7r 4 + 3p 2 



(l+j9sin0 o ) 5 4(l-p 2 ) 5 / 2 ' 



(59) 



hence 



16 7 4 /3 2 



dpp 



{4+p 2y 1 _ {p/f3 y 

(l_ p 2)7/2 



+ 



(4 + 3p 2 ) 



(l-Wy/l-ip/W 



(60) 



which results in 



16 7 4 /3 2 



2/3 2 (2-/3 2 ) 2/3 2 (6 + /3 2 ) 



3(1 -/3 2 ) 2 3(1 — /3 2 ) ; 



16/3 2 



16 7 4 /3 2 3(1 - 13 



2\2 



(61) 
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We perform now the calculation in eq. ( )55|) numerically. Knowing that G(f3, 1) = 1, this 
calculation actually shows the power radiated by iV charges, divided by the power radiated 
by a single charge. The results are shown as function of N for different values of (3 in 
Figure (TUl 




5 10 15 20 25 30 35 



Number of charges N 

FIG. 10: Result of G (eq. (|55p ) as function of the number of charges N, for different values of j3. 
For big N, G goes asymptotically to and as smaller j3 is, G goes faster to 0. 



A. The radiation reaction 

We calculate now the Lorentz force on the charges and the resulting radiation resistance 
power for comparing with the radiated power. We will need the electric field on charge n at 
time t due to charge m at its retarded position at the earlier time t', so that: 
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\v' n (t)-T' m (f)\=c(t-t'). (62) 
Using the expression for r\ in eq. ([1]) we obtain 

~w(t-t) + 2n(n-m)/N~ 



Ad sin 2 



c\t-t')\ (63) 



2 

which is exact for any m and n. By definition, t — t' > 0, so to take the correct square 
root from the left side of this equation, we need to know the connection between m and n. 
To simplify, we restrict: 



0<m<n<N-l, (64) 



for which we obtain 



2rfsin$ nm = c(t-t'), (65) 



where $ nm is defined by 



. ""- f)+ f- mlffl . (66) 
Now we isolate t — t' from eq. fl66|) and set it in eq. fl65l) . obtaining 

$ nm = 7r(n-m)/iV + /3sin$ nm , (67) 

which is an implicit equation, that can be solved by setting a 1st guess $ nm = 7r(n—m)/N 
in the right side of the equation and recalculate $ nm till convergence is obtained. 
Figure ( ITTj) gives the geometrical interpretation of eqs. ()65l)-(l67l). 

Let us look at a simple example of solution for eq. ( 1671) . Say there are 4 charges, so their 
locations at t = are: 0°, 90°, 180° and 270° for charges 0, 1, 2, 3 respectively and let 
us take (3 = 0.7. For calculating the effect of charges 0, 1 and 2 on charge 3, we need to 
calculate $30, $31 and $32, which come out: 2.67259, 2.15479 and 1.48268 respectively, in 
radians. The retarded angle of charge m is 270° — 2$ 3m - see Figure (ITTj) . So translating 
into degrees, we get the retarded angles of —36.256°, 23.079° and 100.097° for charges 0, 1 
and 2 respectively, which are all smaller than the current angles of those charges. 
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Charge m at time t 1 



FIG. 11: Geometrical interpretation for eqs. (|65j) - (|67j ). The charges rotate on the big circle of 
radius d. The retarded distance c(t — t') is the big (red) segment on the 2& mn arc, hence equal to 
2cisin$ mn according to eq. (f65|) . We also see that the big angle 2<fr TO „ (marked in red) equals the 
sum of the angles u(t — t') and 2ir(n — m)/N according to eq. (|66p . Two orthogonal unit vectors 
(green) r(t) and (pit) are drawn near charge n, representing the radial and tangential directions of 
the moving charge. 



It is to be mentioned that the solution $ mn of eq. (167|) is time independent, meaning that 
the angle difference between the current position of charge n and retarded position of charge 
m does not depend on time. Because the charges rotate, the only thing which depends on 
time are the local unit vectors comoving with the charge r(t) and <p(t) - see Figure ( jTTT) . 



0,Q 



ow the electric field on charge n due to charge m at its retarded position is given by 



E, 



q/N 
4vre 



0, 



+ 



cR (1 — 3 ■ R ) 3 



(68) 
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where the 2nd part is the far field which we used in eq. ffTTT) (after replacing a by /3c and 
by c) and the first part is the near field which behaves like 1/R 2 . 

The quantities appearing in eq. (168]) are: f3 m is the retarded velocity of charge m (relative 
to c), R mn is the distance between the retarded position of charge m and the current position 
of charge n (and equals to c{t — t') - see Figure (TTTT) ). and R mn is the unit vector pointing 
from the retarded position of charge m to the current position of charge n. 

We need the field at charge n in local components r(t) and (p{t) (see Figure (TTTTl ) and we 
calculate now all the needed quantities to express the electric field E mn . So we obtain 

R mn - /3 m = r(t) [sin $ nm - /3 sin(2$ nrn )] + Cp{t) [cos $ nm - (5 cos(2$ nm )] (69) 



Rmn x [(R mn - (3 m ) x f3 m ] = w/3(cos $ nm - /3) [r(t) cos $ nm - <^(t) sin $ nm ] (70) 

^ m • R mn = /3 cos $ nm (71) 

and 

Rmn = 2d sin $ nm (72) 

which is easily derived also from Figure ( ITT1) . because i? mn = c(t — £'), which equals to 
2<isin$ nm according to eq. ( 1631) . Putting eqs. (l69l)-(T72l) in eq. ( 1681) we obtain 



F g/^ 1 



4ire Ad? sin $ nm (l - /3 cos $ nm ) 3 

r(t) [sin $ nm - g sin(2$ nm )] + gft) [cos $ nm - g cos(2$ nm )] | 

7 2 sin $„ m 

2/3 2 (cos $ nm - p) [r(t) cos $„ m - ^(t) sin $ nm ]] (73) 

With the aid of this field we will calculate the total force acted on a charge by the other 
charges, but we also need the "self" radiation reaction force. For the most general case, this 



is given by 



3] 
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' self 



{q/Nf 
47ie c 2 



^ 7 2 /3 + 2 7 4 (/3 • P)$ + ■ (3)f3 + 2 7 6 (/3 ■ /3f(3 



(74) 



In our case, of circular motion, the velocity is perpendicular to the acceleration so /3 ■ 
0, we therefore remain with 2 terms 



' self 



(q/Nf 



jj 7 2 /3+^ 7 4 (/3-/3)/3 



(75) 



For the circular motion we know that (3 = —u 2 f3, and using u = v/d = (3c/d, eq. ((75 
reduces to 



(g/AOV 



self 



(76) 



6%e d 2 

showing that the self reaction force is in the direction opposite to the velocity of the 
charge. 

Now we chose n to be the "last" charge, i.e. n = N — 1, hence the restriction in eq. ( )64|) 
holds, and calculate the total force on it, given by the self force plus the force acted by all 
other charges (i.e. the Lorentz force): 



7V-2 



' 7V-1 



' self ' 



(q/N) 



(77) 



m=0 



where vjv-i(t) is the velocity of charge N — 1 and B TOi jv-i is the retarded magnetic field 
on charge N — 1 due to charge m. 

The dumping power on charge iV — 1 is given by Vjv_i(i) • F/v-i, therefore we do not need 
the magnetic part, obtaining 



N-2 

-Pdump 1 charge ^N—li.^) ' FjV— 1 ^N—l 

(78) 

The velocity of the given charge N —1 is in the tangential direction, i.e. \N-iif) = cj3(p(t), 
so only the tangential part of eq. f!73|) affects the power. We obtain 
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-Pdump 1 charge C/3 



(g/AQV 

6ne * d2 m =o 
COS $JV-l,m - /3cos(2$jv_i im 



AT-2 



(,/JV) : 



7 2 sin $ 



JV— l,m 



47re 4ci 2 sin $jV-l,m(l - /9 COS $jV-l, m ) 3 

2/3 2 (cos$jv_i, m - /3) sin^Ar-iJ \ ■ 



(79) 



Because the dumping power on one charge does not depend on time and by symmetry is 
the same for all charges, the total dumping power on the whole system of charges is just N 
times the above: 



P 



dump 



c(3 



<fl A 



N-2 



Qne Nd 2 



r 



1 



m=0 

COS $AT-l,m - COs(2$Af_i 

7 2 sin$Ar_i jm 



A-Ke N Ad 2 sin ® N - hm (l - (3 cos $7v-i,m) 3 



2/3 2 (cos$ 



N—l,m 



(3) sin $ 



N-l,m 



(80) 



The dumping power must be identical with the radiated power, with a minus sign, so to 
compare them, we may factor out — P \ N=1 from eq. ( fTOl ): 



dump 



P\ N=l G ? dump(/^, N), 



(81) 



where 



G ! dump(/3, -^V) 



1 



N-2 



1 



cos $ 7v_i )m — /3 cos(2$ 



N 8JV/3V^ sin^-i,m(l 



N-l,m) 



7 2 sin$jv-i,- 



/3cos $AT_i im ) : 
2/3 2 (cos$ JV _ 1 ,, 



/3) sin $ 



N—l,m 



(82) 



which clearly shows that for N = 1, the sum is 0, so that Gdump(/3, 1) = 1 f° r an Y 
The numerical calculation of eq. ( 1821) shows identical results with those of G calculated from 
eq. ( 1551) . as shown in Figure fTUl 
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IV. ASYMPTOTIC RESULT FOR MANY CHARGES 



To understand how the radiation goes to when the number of charges N goes to infinity, 
one may try to approximate either G from eq. (|55p or Gdump from eq. (|82p . for big N. 
Although Gdump seems more compact, it is more difficult to handle, and we shall develop G 
for large N. 

As mentioned in the previous section (see Figures [3HH]) , for large N the functions Fc and 
Fs tend to be harmonic, hence we may approximate them by the first term of their Fourier 
series. 

For a function a(x) with periodicity X, we specify the Fourier coefficients by 

r x 

A n = a(x) 

JO 

and a(x) is represented by its Fourier series 

1 oo 

a(x) = - ^e l2mx/X (84) 

n=— oo 

For brevity, to refer to the functions Fc and Fs we call them Fc,s (as in eq. ( |34l) ). We 
know those functions have a periodicity of 2ir/N in ip' (see eq. ( 15T1) ). so we define their 
Fourier coefficients: 

Ac,s„ = / Fc,s(^)e- 4iVV V, (85) 
Jo 

and the functions Fc and Fs are expressed as: 

oo 

' n=— oo 

Now we calculate the Fourier coefficients Ac,s n in eq. f )85|) . The integrand is periodic in 
27v/N, so increasing the integration interval to 27r multiplies the result by N, hence we may 
express 

Ac,s n = 1 ^ Fc,s(^)e-^'# / , (87) 
and by using the definitions of Fc,s (eqs. f )20l) and fTST]) ) and the property of 4>k from 
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eq. ( 1501) . we get 

-27T N-l 



Ac,s n = 1 / * ^ fc,s(0 o (^ + 2*klN))e* tw 'dtf/. (88) 

i — n 



fc=0 

We interchange the sum and the integral and change variable tp" = p>' + 2nk/N, obtaining 

N-l <-2w{l+k/N) 
— I J2nk/N 



1 N-l p2w(l+k/N) 

Ac ' s ™ = TrH ***** / fc,s(0 o (^))e-^"V. (89) 



In the above integral, fc,s has a periodicity of 27r and e~ tNnif has a periodicity of 2ir/N , 
therefore the integrand is periodic by 27r. We may therefore move the integration range to 
be between and 2tt, showing that the integral does not depend on k. After renaming ip" 
to tp' we obtain 



1 /N-l \ „2tt p2n 

Ac > s « = a7 ( E e ~ in2Kk ) J kAM<P'))e- iNntp 'd<p' = J o fc,s(0 o (^))e^V, (90) 

because e~ m27Tk = 1 for any k. We see that the n Fourier coefficient of Fc,s is actually the 
Nn Fourier coefficient of fc,s(0o) ; which have a 2% periodicity. This means that if we repre- 
sented the fc,s functions by their Fourier components and evaluated Fc,s = J^^ 1 fc,s(0fc), 
all Fourier components would cancel out except of the nN components, i.e. the 0, N, 2N, 
etc. The Fourier coefficient is 0, because fc,s have DC level (see eqs (I39ti47p ). and clearly 
the second Fourier coefficient of Fc,s, which is the 2N Fourier coefficient of fc,s, is much 
smaller than the first coefficient for large N, as evident also from Figures (J3j)- (IHJ) - 

Therefore, for large N we get the asymptotic Fc,s from its first Fourier coefficient (i.e. 
coefficients 1 and —1), so that we get from eq. ([50]) : 

Fc,s(y/) -> ^ ( Ac, Sl e^ v + Acs.^"^') • (91) 
2tt V / 

Because Fc,s are real, Ac,s_ 1 = Ac,s* and we obtain 

N 

Fc,s(<p') -»■ — lAc^l cos [Np' + arg(Ac,s 1 )] . (92) 

7T 

So we have to calculate the first Fourier coefficients of the Fc and Fs functions, Ac^. 
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From eq. (190]) we get 



Ac, Sl = / fc,s(0o(^'))e-^V, (93) 
Jo 

We change variable to 0o and according to eq. (137|) we have dip'/d(p = — 1 — p sin0 o , so 
we obtain: 

/■(/.o(0)-27r i afAiA 

Ac, Sl = / n f • i , 3 e-^-^-^^0 (-l-p sin0 o ), (94) 

i^o(o) (l+psm0 o ) 3 

where hc,s is an abbreviation for the functions hc(0 o ) and hs(0o), defined in eqs. fl35|) and 
(|36|) . respectively. 

The integrand being periodic on 2tc, we may shift the limits by any value getting 

A c , Sl = /'"J^Le^e^^^o. (95) 
J (l+psm0 o ) J 

We change variable z = exp(z0o) to solve this integral on the complex plane, obtaining: 

Ac, Sl = % I dz . 2 ^[ Z) z N e Nb ^ = -Jdz f(z) (96) 
V Jc \ z + o z - 1) V Jc 

where b is the pure imaginary number defined in eq. fj4~3l . C is the counterclockwise unit 
circle integration contour (see Figure (Q) and the functions lc,s are abbreviations for \c(z) 
and ls(z) defined in eqs. (PfTl) and (1421) . respectively. To ease on further manipulations we 
called the integrand f(z). 

The integrand has two 2nd order poles, 2^ 2, defined in eq. (jUJ) and only z\ lies inside the 
integration contour C - see Figure (Q. In addition there is an essential singularity at z — 0, 
because of the z~ x in the exponent. 

We first calculate the residue at z = zl. Rewriting f(z) 

f(z) = Lc,s(z)z N e N ^ +z - 1 \ (97) 
where Lc,s(z) is defined in eq. f|45|) . we obtain 



Res(/,zl) = hm i/(^ - ^) 2 ] = -f 



lc,s(z) N Nb(z-\-z~ 1 ) 



(98) 



which evaluates to 
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Res(/,zl) 



lc,s'(zi)(zi - z 2 ) - 21c,s'(zi) j^^^-i) ^ lc,s(zi) d 



(zi - z 2 f 



+ 



' z\ — Z2) 2 dz 



gNgNbiz+z- 1 ) 



Z=Z\ 

(99) 

We already showed that the first part is - (see eq. (j460 ). This is because the integral in 
eq. ( |95l) reduces for iV = to the integral in eq. (139]) (up to 2n). So we are left with 



Res(/, zl) 



lc,s(^i) d 
[zi - z 2 ) 2 dz 



lc,s(zi) 
(zi - z 2 ) 



. z N-l e Nb{z 1+ z^) N ^ +b ( Zi _ z -l^ 

(100) 



By using z\ = —l/z 2 and z\ + z 2 = —2i/p = —1/b, we see that this part is too, hence 
the contribution of the pole at z = Z\ is 0. We may therefore exclude this pole from the 
integration contour, and rewrite eq. (|96|) 



2% I , lc,s(3) 
Ac,Sj = — f dz — 

v Jci \ z +0 l z — ly 



V 



n=0 



lc,s(z) 



„Nbz 



(Nb) n 



Ci 



z 2 + b- i z-r 



n\z n ~ N ' 
(101) 

where C\ is the counterclockwise circle of radius smaller than shown in Figure (j9]). 
Inside this integration contour we have only the essential singularity at z — 0, and for 
handling it, we represented the exponent with negative powers of z as a Laurent series. The 
terms n < N are analytic inside Ci, hence contribute to the integral, so by changing the 
summation variable n' = n — (N + 1) and remaining n' to n we obtain 



2i^ (Nb) n+N+1 

Ac > s i ~~ 12 



P 2 ^{n + N + l)\J Cl 



dzLc,s(z)e 



Nbz 



yn+1 



2ij^ (Nb) n+N+1 



dzg(z) 



(102) 

where we used again the definition of Lc,s(z) from eq. (1451) . and the integrand has been 
called g(z), to ease manipulations. We calculate now the residue of g(z) 



Bes(g,0) 



1 d n 
n\ dz 1 



-[Lc,s(z)e 



Nbz] 



,, n\ * — ' \mj \ dz m 

2=0 m= o 



d n ~ n 
dz n ~ 1 



.Nbz 



2=0 

(103) 
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which comes out 



Res(g, 0) = (Nb) n 



{my 



m Jm 



m=0 



'n — m)\m\ dz r 



-Lc,s(z)\ z=0 , 



(104) 



We start with Lc. To handle this derivative we express it as 



d r < 



dz T 



-Lc(z) 



d m+1 



2=0 



dz m+l l _ b -l z _ z 2 



(105) 



2 = 



The last rational function is the Fibonacci polynomials generating function, with argu- 
ment b^ 1 . Hence the result is (m + 1)! multiplied by the m + 1 Fibonacci polynomial. This 
may be directly calculated by factorizing the Fibonacci generating function to obtain 



d" 



dz m 

which results in 



Lc(z) 



d 



m+l 



2=0 



Vb- 2 + 4 dz m+1 



1 — z 1 1 z 1 — z 2 1 z 



2=0 



(106) 



dz r 



-Lc(z) 



(m + l) 



_1\ m+l 



2=0 



Vb~ 2 + 4 



or explicitly 



d r ' 



dz r 



+ y^TI)^ 1 - (fe- 1 - y^+4) 



= (m + 1) 



2=0 



+ 4 



The Ls case is handled similarly. We express: 



dT 



dz T 



= ib 



d 



m+l 



2 - b~ l z 



2 = 



dz m+l l _ b -l z _ z 2 



2 = 



m+l 



(107) 



(108) 



(109) 



The last rational function is the Lucas polynomials generating function, with argument 
b^ 1 . Hence the result is (m + l)! multiplied by the m + l Lucas polynomial. This may be 
directly calculated by factorizing the Lucas generating function to obtain 



dz r 



: Ls(z) 



ib 



d 



m+l 



2 = 



dz m+1 



1 — z l l z ~^ 1 — z 2 1 z 



(110) 



2=0 



which results in 



dz r 



2 = 



ib(m + l)\ (z;T +1 + fe" 1 ) 



-nm+l 



(111) 
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or explicitly 



2»n+l 



ib(m + 1)! 



z=0 



m+l 



(112) 



By using eqs. f II 8 [) . f ll04|) and (11021) and replacing 6 = — 0.5ip we obtain a closed form 
expression for Aci 



Aci 



iJV+l 



-l) n (A/2) 



2n+N+l 



n=0 



n + A + 1)! 

m+l 



E 

m=0 
1 - 



(-l) m (m + l) p 2(n-m)+JV-l 

(A/2) m (n-m)! 

m+l~ 



l—p x 



(113) 



and for Asi 



As = -27r(-i) N +H V (-l) n (^/2) 2 " +JV+1 ^ (-l) m (m + 1^(— )+iV-i 
" ^ (n + N+l)\ ^ 



n=0 



m=0 



m+l 



(A/2) m (n -m)! 

m+l 



1 — 



;ii4) 



We remark that if A is a multiple of 4, the angle of Aci is 90°, and each increment of N 
removes 90° from the angle of Aci. Also we see that the angle of Asi is always bigger by 90° 
than the angle of Aci and this is visible in Figures ([3]), fl5]) and ©. We shall name those 
angles ip c and ip s in the following calculations. 

So by setting eqs. (192]) in ( 1551) . and inverting the integration order between dp and dip' 
we obtain: 



G((3, N -»■ oo) 



47tA 7 4 /3 2 



o V 1 " (P// 3 ) 2 Jo 



2n/N 



dp' cos 2 (Np' + p s ) + 



d PP i-\A Cl \j vi-(p//3) 2 y o ■ 



d<^' cos ( A<£>' + <£> c 



(115) 
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The dip' integrals result in half the integration interval, i.e. 7r/N, so after simplifying we 
obtain 



G(J3,N -> oo) 



N 



dpp 



|AsrP 



VI " (p/0) : 



+ |A Ci |V1-(p//3) 2 • ( 116 ) 



Figure f lT2|) shows the asymptotic results of G for large N according to eqs. fll!6j) . fll 13j) 
and (11141) . compared with the exact result from eq. (|55|) . 
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FIG. 12: Asymptotic results for G according to eq. (I116P versus exact results according to eq. (|55l) 
as function of the number of charges N, for different values of /3. The asymptotic approximation 
is accurate even for a few number of charges. As /3 gets bigger, the difference between asymptotic 
and exact is more visible for small N. 



The asymptotic result is much easier calculable than the exact one, and does not require 
to solve for each step the implicit equation (13"2"j) . but is still not given by a simple formula. 

We will calculate in the next subsection an asymptotic expression for small (3, i.e. G((3 — > 
0, — » oo), and for this case one arrives to a simple formula, as we shall see below. 
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A. Asymptotic result for many charges and low velocity 

For small /3, |p| <C 1, hence -^/l — p 2 in the denominator of eq. (11131) may be set to 1, 
and 1 + yl — p 2 ~ 2, neglecting 1 — yl — p 2 in eqs. (11131) and (11141) . So we obtain 



Ac 1 | w _ MJ = -27r(-») 2. (n + A + 1)! ^ (JV/2)™(n -m)! 2 (H7) 

n=0 K ' m=0 \ I J \ ) 

and 



Asilip^o = *ACi||p|_^ - ( 118 ) 

Rearranging eq. (I117P we obtain 



Aci| 



-4tt( z; p 2^ r n + A r + i)i ZWr,.- 



n=0 



m=0 



(n — m)\ \Np 



(119) 



and defining K = —4/ (Ap 2 ), one may express the sum over m in eq. fl 1 1 9 [) as a derivative 
with respect to K as follows 



A 



m+l 



m=0 



(n — m)! 



and by changing the summation variable m' = n — m, this is written as 



j^n—m' + l 

9k ~< — = 9k 



Rn+1 A (1/A)™' 



m'=0 



7ft ! 



We may sum exactly the last sum over m' to obtain 



(120) 



121) 



2^ -/i = ex P(VA) — Z1 — * — - « exp(l/A), 



m'=0 



m 



(122) 



where T with 2 arguments is the incomplete gamma function. We are interested in small 
|p|, so for |1/A| = Ap 2 /4 <C 1, we obtain the approximated result given in eq. (I122p . which 
means that for small argument, the exponential series in eq. f II 2 2 [) needs very few terms to 
converge to an exponent. So continuing the calculation started in eq. (I12ip we obtain 
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d K [K n+1 ex V (l/K)}] = K n ex V (l/K)[n+l-l/K] « (n + l)K n = (n + l) , (123) 

where the last approximation used again the fact that \1/K\ <C 1. Now using the result 
from eq. f 11 2 3 j) in eq. (I119p . we obtain 

oo 

A Cl | bM = -4 7 r(-i)^- 1 (iV/2)^ 1 £ (n + N + l)\ Nn - (124) 

n=0 ^ ' 

The above may be summed exactly, obtaining 



^ n + l _ N N+3 T(N) - e N (T(N + 2)T(N + 1,N) - NT(N)T(N + 2,N)) 

^ (n + N + 1)! ~ jV( 2 +*0r(iV)r(iV + 2) ' 

(125) 

For large N,T(N + 1,N) « ±r(iV + l) and T(iV + 2, iV) « |r(JV+2), hence the expression 
multiplying the exponent in eq. fl 1 2 5 [) tents to 0, remaining with 



oo 

n 



y N « ~ N ' [2+N) i NN+3 nN)} = N _ 1_ 

^(n + N + l)\ ~ T(N)T(N + 2) T(N + 2)~N\' V ; 



n=0 

Putting eq. f TT26l) in <II21)> we obtain 



-i 1 



A Cl | bh0 = -47r(-i) iV +W(A r /2) iV+1 ^ « (-z) Ar+1 v / 27iV(e/2) V" 1 , (127) 

where the last expression has been obtained by using the Stirling approximation, and 
Asi | |p|— »-o = iAci | ipi— >-o according to eq. (11181) . 

Now we may perform the integral in eq. (11161) . which becomes for small (3 



4vrV/3 2 Jo {y^l-ip/py 

(128) 

Because — > 0, we may set 7 = 1, and we obtain 
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(129) 



For large N, T(N + 3/2)/T(N) « iV 3/2 , so we obtain 



G(f3 ->■ 0,iV ->■ oo) 



3N(e/2) 2N P 2N (N + l)yft 3(e0/2) 



2iV 



AT 



(130) 



2vr/3 2 2iV 3 / 2 ' 4v^F/3 2 

Figure [13] shows the asymptotic results calculated with eq. ( 1130|) versus exact results 
according to eq. (1551) . 
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FIG. 13: Asymptotic results for G according to eq. (|129p versus exact results according to eq. ([55]) 
as function of the number of charges N, for different values of (3. For (3 = 0.05 the asymptotic 
result is completely indistinguishable from the exact result, while for (3 = 0.1 and f3 = 0.15 they 
are almost indistinguishable. For (3 > 0.2 the asymptotic approximation becomes inaccurate. 



It is to be mentioned that the case of small (3 and large iV analyzed here fits the situations 
of currents in conducting materials or ion drift currents, mentioned in the introduction. In 
a conducting loop, the number of charges may be of order of 10 23 , and (3 may be of order 
10~ 12 , so that f3 2N results in completely unmeasurable radiated power. In a ion drift device, 
the number of charges may be of order of 10 10 and (3 may be of order 10 -6 , so that the 



32 



radiated power is somehow bigger than for the conducting loop, but still unmeasurable. 

V. CONCLUSIONS 

The purpose of this work was to learn how the radiation from discreet charges vanishes in 
the continuum steady state limit. We used a canonical configuration of charges in uniform 
circular motion, uniformly spread around a circle. 

We found that the log of the power decreases almost linearly with the increase in the 
number of charges, and arrived to a close form solution to calculate this power if the number 
of charges is big - see Figure (fl~2l) . 

Specifically, for low velocities, we derived a simple expression for the radiated power. It 
shows that the radiated power is governed by ft at the power of twice the number of charges 
- see eq. ( I130p and Figure (ITBT) . explaining why the radiation in all the cases considered as 
DC is unmeasurable. 
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